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Domestication of fruit and vegetables resulted in a huge diversity of shapes and sizes of 
the produce. Selections that took place over thousands of years of alleles that increased 
fruit weight and altered shape for specific culinary uses provide a wealth of resources 
to study the molecular bases of this diversity. Tomato (Solanum lycopersicum) evolved 
from a wild ancestor (S. pimpinellifolium) bearing small and round edible fruit. Molecular 
genetic studies led to the identification of two genes selected for fruit weight: FW2.2 
encoding a member of the Cell Number Regulator family; and FW3.2 encoding a P450 
enzyme and the ortholog of KLUH. Four genes were identified that were selected for fruit 
shape: SUN encoding a member of the IQD family of calmodulin-binding proteins leading 
to fruit elongation; OVATE encoding a member of the OVATE family proteins involved in 
transcriptional repression leading to fruit elongation; LC encoding most likely the ortholog 
of WUSCHEL controlling meristem size and locule number; FAS encoding a member in the 
YABBY family controlling locule number leading to flat or oxheart shape. For this article, 
we will provide an overview of the putative function of the known genes, when during 
floral and fruit development they are hypothesized to act and their potential importance in 
regulating morphological diversity in other fruit and vegetable crops. 
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INTRODUCTION 

Angiosperm plants vary tremendously in morphological traits 
related to their reproduction. The floral appearance is driven 
by evolutionary aspects of the pollination syndrome whereas 
distinct dispersal modes drive the evolution of phenotypes asso- 
ciated with the fruit. In natural settings, the main functions 
of the fruit are to protect the developing seeds and to act as 
a dispersal agent. The onset of the change to an agricultural 
lifestyle, approximately 10,000 years ago, provided strong selec- 
tion pressures on the fruit of incipient vegetable and fruit crops. 
The selections made by early farmers offer a great opportu- 
nity to identify the molecular basis of a range of phenotypic 
traits, especially those related to fruit morphology and flavor. 
For example, selections against bitter taste resulted in palatable 
eggplant and cucumber (Wang et al., 2008; Qi et al., 2013). 
Yet, the underlying principle for nearly all cultivated vegetable 
and fruit crops was the selection for larger and more nutri- 
tious fruits featuring a variety of shapes (Paran and Van Der 
Knaap, 2007; Pickersgill, 2007; Meyer and Purugganan, 2013) 
(Figures 1A-C). The larger fruit became more nutritious as a 
result of the increase in the edible and fleshy part of the fruit 



at the expense of the seed part for most domesticated fruits and 
vegetables. 

The focus of the "hypothesis and theory" article is to summa- 
rize the current knowledge on the function of genes that change 
tomato fruit weight and shape resulting from domestication and 
diversification process. The focus on tomato is based on the exten- 
sive research that resulted in the cloning of six fruit shape and 
weight genes from this species in recent years. The predicted func- 
tion of these genes will be discussed in the context of the phases of 
development where we hypothesize the impact of the mutant alle- 
les is most critical. It is important to recognize that the mutations 
are not often resulting in complete nulls, i.e., a loss-of-function 
allele. Thus, the complete repertoire of functions of the tomato 
fruit shape and weight genes may not be apparent from the phe- 
notype observed in the natural mutants. We will propose the 
pathways in which the shape and weight proteins function. We 
will also include the molecular basis of the underlying mutations 
that gave rise to the derived alleles and demonstrate that inver- 
sions, duplications, as well as single nucleotide polymorphisms 
(SNPs) in promoters and coding regions underlie the phenotypic 
diversity of the tomato fruit. 
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FIGURE 1 | Diversity in tomato fruit shapes. (A) Tomato varieties carrying 
the ovate mutation result in obovoid and ellipsoid fruit shapes. Size bar = 
3 cm. (B) Tomato varieties carrying the sun and Ic mutation result in a long 
fruit shape. The oxheart shaped tomato also carries fas in addition to sun 
and Ic. Size bar = 3 cm. (C) Tomato varieties carrying fas and Ic result in a 
flat fruit with many locules. The wild type represents the fruit from an 
ancestor of cultivated tomato, S. pimpinellifolium LA1 589. Size bar = 3 cm. 
(D) The axes of growth: proximal (closest to the stem) to distal (farthest 
away from the stem); medio-lateral; and the adaxial (closest to the 
meristem) to abaxial (farthest from the meristem). The fruit tissues that 
give structure to the organ are highlighted as the pericarp, septum, 
columella and locule. Not indicated is the placenta which is the tissue 
extending from the columella and surrounds the seeds. 



OVERVIEW OF TOMATO DEVELOPMENT 

Even though the fruit is a terminal structure that forms relatively 
late in the plant's lifecycle, the formation of this organ and the 
parameters that determine its final dimensions are rooted much 
earlier in the plant's lifespan. Therefore, it is important to view 



tomato fruit development in the context of overall plant devel- 
opment starting after germination. Plant growth in tomato and 
other Solanaceous plants is characterized by a sympodial shoot 
architecture where after formation of 8-10 leaves, the shoot api- 
cal meristem (SAM) terminates into the inflorescence meristem 
(IM), and growth continues from lateral meristems called sym- 
podial meristems (SYM). Meanwhile, the IM terminates into 
the floral meristem (FM) generating the flower (Schmitz and 
Theres, 1999). The tomato inflorescence also features a sympo- 
dial structure since a new IM emerges simultaneously from the 
flank of the first FM, terminating again in the second FM on 
the inflorescence and so on (Figure 2A). This growth pattern is 
referred to as cymose and results in a zigzag of flowers on a 
tomato inflorescence (Welty et al., 2007; Lippman et al., 2008; 
Castel et al, 2010). In most angiosperm species, FMs give rise 
to four whorls: the sepals, petals, stamens and carpels. Organ 
identity genes play critical roles to ensure that carpel primor- 
dia arise from specified founder cells within the FM (Causier 
et al, 2010) (Figures 2A,B). In addition to cell specification, the 
establishment of the boundaries between and within the pri- 
mordia is required to ensure that the appropriate identities and 
division patterns are initiated and maintained throughout gynoe- 
cium growth (Dinneny and Yanofsky, 2005; Balanza et al., 2006; 
Girin et al., 2009). This step is critical to lay the foundation of 
growth of the organs along three axes: the proximal-distal, the 
medio-lateral and the abaxial-adaxial axis (Figure ID). A mature 
tomato gynoecium coincides with flower opening which marks 
the anthesis and pollen release stage (Xiao et al, 2009). Following 
pollination and fertilization of the ovules, fruit development is 
initiated which is marked by a rapid increase in cell prolifera- 
tion followed by cell enlargement (Gillaspy et al., 1993; Xiao et al, 
2009) (Figures 2G,H). In most fruit tissues such as the pericarp, 
cell division ceases 5-10 days after anthesis and growth of the 
fruit continues by extensive cell enlargements that last for three 
to 5 weeks until the fruit ripening stage (Gillaspy et al., 1993; Xiao 
et al, 2009). 

CRITICAL DEVELOPMENTAL STAGES BEFORE ANTHESIS AT 
WHICH THE FINAL SHAPE AND WEIGHT OF FRUIT IS 
REGULATED 

The final dimensions of the fruit are regulated during multiple 
stages throughout the development of the plant. These stages 
occur before and after anthesis, and may be initiated as early 
as in the SAM. Thus, the first stage of regulation of the final 
fruit dimensions is likely to occur in the meristems as a result 
of their size (Figure 2A; Table 1). Since the gynoecium is a ter- 
minal structure, the size of the FM may impact the number of 
cells that are specified to form a carpel primordium as well as 
the number of primordia (Szymkowiak and Sussex, 1992; Clark 
et al, 1993; Taguchi-Shiobara et al, 2001; Suzaki et al., 2004). Cell 
identity and the positioning of organ primordia per se however 
are not controlled by the size of the FM. Therefore, the second 
stage of regulation is likely controlled by the organization within 
the meristem which relates to where and how often in the meris- 
tem the cells that are destined to become carpel primordia arise 
(Figures 2A,B). Similarly as for leaf primordia initiation, local- 
ized auxin maxima controlled by the auxin efflux protein PIN1 
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FIGURE 2 | Tomato floral and fruit development in S. pimpinellifolium 
LA1589. (A) The IM terminates into the FM, which develops into a floral 
bud 2 dpi. The SYM is emerging from the flank of the youngest leaf 
primordium. (B) Emergence of carpel primordia 5 dpi. (C) Growth of 
carpels and formation of two locules 6 dpi. (D) Growth of the style. (E) 
Formation of the ovules 10dpi. (F) A transverse section of the pericarp 



at anthesis. (G) Transverse section of the pericarp 5dpa. (H) Transverse 
section of the pericarp 10dpa. SYM, sympodial meristem; IM, 
inflorescence meristem; FM, floral meristem; dpi, days post floral 
initiation; dpa, days post anthesis; p, petal; st, stamen; ca, carpel; sy, 
style; ov, ovule; exo, exocarp; endo, endocarp; meso, mesocarp. Size bar 
is 50p,m. 



Table 1 | Developmental phases proposed to control fruit shape and weight. 

Critical regulatory Developmental event 1 Landmark 1 Cellular events in the ovary or fruit Days after Stage-specific 

phases of fruit meristem fruit shape and 

shape and weight initiation weight genes 



Phase 1 

Phase 2 
Phase 3 
Phase 4 
Phase 5 

Phase 6 

Phase 7 



Inflorescence and floral meristem 
formation 

Floral meristem organization 
Gynoecium initiation 
Gynoecium growth 
Anthesis 

Fertilization and 4-16 cell stage 
embryo 

Globular to coiled stage embryo 



Floral landmark 1 

Floral landmark 1 
Floral landmark 5 
Floral landmark 6-9 
Floral landmark 10 and Flower opening 
fruit landmark 1 

Fruit landmark 2-3 Cell proliferation 



Cell number, size of the stem cell niche 

Cell identity and boundary information 
Cell proliferation and enlargement 
Rediffentiation of tissue types 



Fruit landmark 4-7 



Cell enlargement 



1 

5-6 
8-16 
19 

20-25 

25-39 



LC/FAS/CNR 

FAS/CNR 
OVATE/SUN 



SUN/KLUH 



1 From Xiao et al. (2009). 



and the expression of PLETHORA/ AINTEGUMENTA transcrip- 
tion factors are thought to control floral organ positioning 
(Benkova et al, 2003; Krizek, 2011; Van Mourik et al, 2012; 
Hofhuis et al, 2013). The areas of low auxin coincide with the 
boundaries between primordia which are also tightly controlled 
processes (Nahar et al., 2012; Zadnikova and Simon, 2014). 
Misalignment during this stage would result in changes in final 
fruit morphology. The third stage is the phase that transmits 
positioning information to gynoecium growth (Figures 2B,C, 
Table 1). During this phase, the three axes of growth have 
been specified along which cell proliferation and enlargement 
occurs (Dinneny et al., 2005; Ostergaard, 2009). Cell prolifer- 
ation, which is characteristic of this stage, consists of the rate 



and duration of the cell divisions within the developing ovary 
impacting final organ dimensions (Figure 2C). Also critical are 
the differential rates and duration of cell division within dis- 
tinct tissues in the developing ovary, resulting in alternatively 
shaped fruit. For example, ovary and fruit length is determined 
by the degree of growth along the proximal-distal axis whereas 
width is determined by the degree of growth in the medio- 
lateral axis (Figure ID). The degree of the pericarp thickness 
and other internal tissues is determined along the abaxial-adaxial 
axis. Therefore, enhanced cell divisions preferentially along one 
axis of growth are proposed to lead to a different shape fruit as 
opposed to enhanced cell divisions along all three axes of growth. 
The fourth stage occurs concomitantly with the third stage which 
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is the continued specification of new tissue types through reac- 
tivation of the meristematic potential leading to the formation 
of many tissue types (Girin et al., 2009) (Figures 2D,E). Along 
the proximal-distal axis, the gynoecium develops two additional 
regions: the stigma and style. Along the medio-lateral axis, the 
ovary develops the placenta, ovules and transmitting track tissues. 
Along the abaxial-adaxial axis the ovary continues to maintain 
the polarities within the different tissues such as the pericarp, 
septum, placenta and ovules. The reinforcement to maintain the 
different zones is mediated by transcription factors in conjunc- 
tion with boundary genes (Heisler et al, 2001; Nahar et al., 
2012). 

CRITICAL DEVELOPMENTAL STAGES AFTER ANTHESIS AT 
WHICH THE FINAL SHAPE AND WEIGHT OF FRUIT IS 
REGULATED 

The anthesis/pollination/fertilization phase marks the end of 
ovary development and the beginning of fruit development. 
Lack of or poor fertilization leads to changes in fruit shape and 
reduced weight, marking the fifth phase. Aborted fruit is ter- 
minal and should not been considered to be part of phase 5. 
The first stage post-anthesis is the sixth phase proposed to cor- 
respond to the cell proliferation stage, a rapid increase in cell 
division throughout the developing fruit that follows immedi- 
ately after fertilization (Gillaspy et al, 1993; Xiao et al, 2009) 
(Figures 2F,G). As in the ovary, this stage is comprised of dif- 
fering cell division rates and duration in the tissues of the fruit 
that would greatly impact final fruit shape. The seventh and final 
stage is proposed to be cell enlargement which impacts overall 



fruit size the most (Figure 2H). Cell enlargement is regulated dif- 
ferentially in the various tissues within the fruit, and rates and 
duration determine the final fruit dimensions. For example, the 
columella and placenta tissues contain more large cells than the 
pericarp. Additionally within the pericarp, the exocarp cells (con- 
stituting the epidermis) are very small whereas the mesocarp cells 
are large (Figure 2H). 

TOMATO FRUIT WEIGHT AND SHAPE ALLELES ACTING 

PRE-ANTHESIS 

L0CULE NUMBER 

LOCULE NUMBER (LC) controls the number of carpel primor- 
dia and a mutation results in a fruit with more than the typical 
two to three locules (Barrero et al., 2006; Munos et al, 2011). 
Increases in locule number often lead to a flat fruit of a larger 
size and the mutation is common in beefsteak tomato and toma- 
toes on the vine (Munos et al., 2011; Rodriguez et al, 2011) 
(Figure 1C). Since carpel primordia arise early in floral develop- 
ment, it is likely this gene functions in regulating meristem size 
and/or in the initiation of organ primordia. The locus was fine 
mapped to a 1608 bp region located between a putative ortholog 
of WUSCHEL (WUS) (annotated gene ID Solyc02g083950, avail- 
able at http://solgenomics.net/) and a WD40 motif containing 
protein (Solyc02g083940). Further association mapping led to the 
identification of two single nucleotide polymorphisms located 
1080 bp downstream of the putative tomato ortholog of WUS 
(Munos et al., 2011) (Figure 3). WUS encodes a homeodomain 
transcription factor that is required for maintaining the stem cell 
identity in the SAM (Mayer et al., 1998; Clark, 2001). The WD40 
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FIGURE 3 | The molecular basis of tomato fruit shape and weight 
variation. (A) Genome structure of the fruit shape and weight loci and the 
underlying mutations. Red box indicates the coding region of a functional 
gene whose regulation is altered by the mutation (denoted by X). Pink 
indicates a loss-of-function mutation of the gene. The size of the loci are not 
drawn to scale. (B) Protein features of the fruit shape and weight proteins. 
The box represents the coding region. The most important domains are listed 



as red boxes. IQ67, CaM binding domain of 67 amino acid and containing IQ; 
OFP Ovate Family Protein motif of unknown function; HLH, YABBY type of 
DNA binding domain featuring a helix-loop-helix structure; HD, DNA binding 
homeodomain of the helix-loop-helix-turn-helix structure; WUS, essential for 
proper functioning of WUSCHEL; EAR, transcriptional repressor function; 
PLAC8, similarity to the placenta-specific gene 8 protein; CYP450, 
cytochrome P450. Size bar = 50 amino acids. 
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containing motif protein belongs to a large family involving 
in diverse functions ranging from signal transduction to tran- 
scriptional regulation (Ullah et al., 2008). Increased expression 
of WUS in Arabidopsis leads to increased floral organ num- 
ber, which is similar to the phenotype found in the Ic mutant 
(Mayer et al, 1998; Clark, 2001). Therefore, based on the pre- 
dicted function SIWUS is the most likely candidate to underlie 
Ic, impacting the first phase that regulates the final dimen- 
sion of the tomato fruit (Figure 2A and Tables 1, 2). Similar 
to Arabidopsis, SIWUS is expressed in the youngest floral buds 
and the shoot apex and virtually undetectable in other tomato 
tissues (Figure 4A). Its expression is also high in the IM/FM tis- 
sues, decreasing very rapidly as floral development progresses 
(Figure 4B). 

WUS is critical in the regulation of the stem cell popula- 
tion size in all meristems, yet the Ic mutation itself does not 
lead to dramatic changes in SIWUS gene expression compared 
to wild type (Munos et al, 2011). Therefore, the high locule 
number phenotype is likely due to subtle changes in expres- 
sion that were not captured by the method of gene expression 
quantification. WUS positively regulates the expression of the 
MADS box transcription factor AGAMOUS (AG) (Lenhard et al., 
2001; Lohmann et al., 2001) and AG is critical in determining 



stamen and gynoecium identity (Yanofsky et al, 1990). Therefore, 
WUS-induced expression of AG links meristem activities to organ 
identity processes. AG in turn down-regulates expression of WUS 
providing the mechanism for changing stem cell identity of the 
remaining FM to carpel identity (Lohmann et al, 2001; Liu et al., 
201 1). In Arabidopsis, WUS down-regulation is mediated by two 
downstream CArG cis-regulatory elements to which AG binds, 
resulting in the epigenetic silencing of WUS (Tilly et al., 1998; 
Liu et al, 2011). Intriguingly, the two SNPs located downstream 
of tomato WUSCHEL are located in a putative tomato CArG cis- 
regulatory element (Figure 5A). This suggests that the Ic mutation 
causes a loss-of-function regulatory element permitting higher 
expression of SIWUS and maintenance of a larger stem cell pop- 
ulation resulting in increased locule numbers. Furthermore, this 
finding implies that the Ic mutation acts at the transition from 
stem cell identity to carpel identity acting just prior to the stage 
shown in Figure 2B. Other critical components of the WUS sig- 
naling pathway are provided by the CLAVATA (CLV) proteins 
(Clark, 2001; Brand et al, 2002; Lenhard and Laux, 2003). In par- 
ticular, the WUS and CLV3 feedback loop is tightly linked to the 
regulation of meristem size in Arabidopsis (Schoof et al., 2000), 
suggesting that members of the CLV pathway may be involved 
in the regulation of tomato meristem size and its organization 



Table 2 | List of genes controlling fruit weight and shape variation in tomato. 



Locus/QTL 


Underlying 
gene ID 


Putative cellular/molecular 
function and length of the protein 


Timing of the impact 
on morphology 


Most likely cause of 
allelic variation 


References 


fw2.2 


Cell number 
regulator (CNR) 
Solyc02g090730 


Increased expression is associated 
with reduced cell division. May permit 
transport across membranes. Protein 
may be located at the 
plasmamembrane and contains a 
PLAC8 domain including two putative 
transmembrane motifs. 163aa 


Phase 1 or 2 
(Figures 2A,B) 


SNP in the promoter 
of the gene 


Frary et al., 2000; Guo 
eta!., 2010 


fw3.2 


KLUH 

Solyc03g114940 


A cytochrome P450 of the 78A class 
and the likely ortholog of AtKLUH. 
Hypothesized to synthesize a mobile 
signal. Substrate unknown. 516aa 


Phase 5, (Figure 2G) 


SNP in the promoter 
of the gene 


Anastasiou et al., 2007; 
Chakrabarti etal., 2013 


Ic 


WUSCHEL 
Solyc02g083950 


Homeobox domain protein. Required 
to maintain stem cell identity in 
meristems. 73 aa 


Phase 1, (Figure 2A) 


Two SNP located 
downstream of 
WUSCHEL 


Mayer etal., 1998; 
Munos et al., 2011 


fasciated 


YABBY2 

Solyc11g071810 


Transcription factor involved in organ 
polarity and meristem organization. 
177 aa 


Phase 1 or 2, 
(Figures 2A,B) 


Gene knock out by a 
294 kb inversion with a 
breakpoint in the first 
intron of YAB2 


Cong etal., 2008; Huang 
and Van Der Knaap, 
2011; Huang etal., 2013 


ovate 


OVATE 

Solyc02g085500 


Increased expression is associated 
with shorter plants and plant organs. 
May be a repressor of transcription. 
Contains the OFP domain. 352 aa 


Phase 3, (Figures 2B,C) 


Premature stop codon 
in an exon associated 
with a mutant 
phenotype 


Liu et al., 2002; 
Hackbusch et al., 2005; 
Huang etal., 2013 


sun 


SUN 

Solyc10g079240 


Increased expression is associated 
with elongated fruit. Positive regulator 
of growth. Contains the IQ67 motif 
that binds calmodulin. 421 aa 


Phase 3 and 6, 
(Figures 2B,C,G) 


Interchromosomal 
gene duplication 
mediated by the 
transposon Rider 


Abel et al., 2005; Xiao 
et al., 2008; Huang 
etal., 2013 
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FIGURE 4 | Expression analysis of the six fruit shape and weight genes 
in different tissues and at different developmental stages. Samples were 
collected from the S. pimpinellifolium accession LA1589. The data was 
obtained from 3 to 4 biological replicate RNA samples that were sequenced 
using the HiSeq2000 lllumina sequencing technology (Huang et al., 2013) 
(http://ted.bti.cornell.edu/cgi-bin/TFGD/digital/home.cgi). The expression was 
normalized using the reads per kilobase per million mapped reads of each 
gene model (RPKM). (A) The root, hypocotyl (hypo), cotyledon (cotyl) and 
shoot apex including the SAM (apex) were collected from the same seedlings 



germinated in petri dishes. All other tissues were collected from mature 
plants grown in the greenhouse (Huang et al., 2013). yl, young leaves; ml, 
mature leaves; yfb, young floral buds from 10 dpi and younger; ant, whole 
flower at anthesis; 10 and 20 dpa, developing fruit 10-20 days after anthesis; 
break, breaker stage fruit which is immediately before turning color. (B) 
IM/FM, 2, 4, 6 dpi flower buds that were fixed in RNAIater solution. The 
tissues were hand-dissected using a dissecting scope prior to RNA isolation. 
Each replicate out of 3 is represented by 100-150 samples that were pooled 
prior to RNA extraction. 
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FIGURE 5 | The effect of Ic and fas loci on locule number in tomato. 

(A) Alignment of the wild type (accession JF284938) and mutant 
(JF284939) LC allele sequences with the canonical MADS box 
transcription factor CArG1 binding sequence (Tilly et al., 1998). The two 
mutations in LC reduce the alignment to the consensus sequence. (B) 



The effect on locule number in the fas, Ic and the double fas/lc NIL 
compared to wild type (WT). The number below the NIL indicates the 
average locule number from over 40 fruit evaluated each from 5 plants. 
The increase in locule number in the double NIL indicates synergistic 
interactions of the two mutations. Size bar = 1 cm. 



leading to changes in locule number and the final shape of the 
fruit. 

FASCIATED 

The mutation in FASCIATED (f or fas) leads to increases in 
locule number with more pronounced effects on locule number 



than Ic (Lippman and Tanksley, 2001). fas is found in certain 
heirloom tomatoes and a few commercially grown beefsteak 
varieties (Rodriguez et al., 2011) (Figures 1B,C). In addition to 
increased locule number, the fas mutation results in increased 
number of all floral organs (Lippman and Tanksley, 2001; Barrero 
and Tanksley, 2004). Significant epistatic interactions have been 
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detected between Ic and fas (Lippman and Tanksley, 2001; Barrero 
and Tanksley, 2004), suggesting that both genes act together by 
co-regulating a core pathway that controls locule number. FAS 
was fine mapped to the bottom of chromosome 11 and, con- 
trary to previously reported results, the mutation resulted from a 
294 kb inversion with one of the breakpoints in the first intron of 
a member of the YABBY family creating a null mutation (Huang 
and Van Der Knaap, 2011). This YABBY member, SIYABBY2 
(YAB2) is considered to underlie fas (Cong et al., 2008) (Figure 3). 
Compared to any other fruit shape or weight gene, YAB2 expres- 
sion is very high in cotyledons, shoot apex, young leaves, young 
floral buds, and anthesis stage flowers (Figure 4A). In IM/FM 
and developing floral buds, its expression is relatively low in 
the meristem but increases in flower buds 6 days after initiation 
(Figure 4B). 

The YABBY family of transcription factors is known to con- 
trol the abaxial-adaxial polarity of SAM, IM, and FM, while also 
specifying the cell fate of the abaxial region in lateral organs. 
YABBY proteins function redundantly with other polarity pro- 
teins and are required to establish the proper boundaries within 
the meristem and developing organ primordia (Bowman and 
Smyth, 1999; Bowman et al, 2002). Moreover, YABBYs have been 
shown to impact the signaling from lateral organs to the meris- 
tem and coordinately maintain the normal growth of meristem 
in Arabidopsis and rice (Goldshmidt et al, 2008; Tanaka et al., 
2012). Because of the function of YABBY family proteins and its 
expression pattern, we consider that FAS is controlling the sec- 
ond stage of final fruit size and shape regulation by impacting 
meristem organization and boundary information (Figure 2B, 
Table 2). However, because of its epistatic interaction with LC, it 
is also possible that FAS impacts meristem size as well as organiza- 
tion (Figure 2A). The details of how YAB2 impacts locule number 
are not well understood. 

Of the two loci controlling locule number, Ic and fas, the for- 
mer mutation is much more widespread in the tomato germplasm 
than the latter while the latter has a more dramatic effect on locule 
number resulting in up to countless locules per fruit (Munos 
et al, 2011) (Figures 1B,C). In near-isogenic lines (NILs) using 
the wild species LA1589 as the background, the impact of these 
two genes on locule number is much less dramatic (Figure 5B), 
supporting the notion that in the cultivated background mod- 
ifiers of these mutations exist. Further genetic analyses would 
reveal the molecular nature of those modifiers. The epistatic 
interaction between the two loci is clearly evident in the wild 
species background as locule number increase in the double NIL 
is higher than the sum of locule number found in the single NILs 
(Figure 5B). 

OVATE 

The shape of many ellipsoid and obovoid varieties such as those 
found in grape tomato is controlled by the gene that regu- 
lates fruit elongation, OVATE (Ku et al., 1999; Liu et al., 2002; 
Rodriguez et al., 2011) (Figure 1A). The gene was fine mapped 
to chromosome 2 and the mutation resulted in a premature stop 
codon in a newly defined class of plant proteins, Ovate Family 
Proteins (OFP) (Liu et al., 2002; Hackbusch et al, 2005) (Figure 3 
and Table 2). The expression of wild type OVATE is the highest 



in the shoot apex, youngest floral buds and breaker stage fruit 
(Figure 4A). Additionally, even though OVATE expression is the 
highest in the IM/FM, expression is reduced by only ~30% in 
flower buds 2, 4, and 6 days after initiation (Figure 4B); the latter 
stage corresponds to the stage shown in Figure 2C. Not all tomato 
varieties that carry the ovate mutation display an elongated shape 
which led to the mapping of two suppressor loci, sovl and sov2, 
on chromosomes 10 and 11, respectively (Rodriguez et al., 2013). 
These suppressors are thought to play important roles in the reg- 
ulation of shape mediated by the OVATE pathway. OVATE does 
neither affect floral organ identity, FM organization nor floral 
organ number (Liu et al., 2002). Instead, OVATE appears to have 
a specific role in the regulation of anisotropic growth along the 
proximal-distal axis at the proximal end of the fruit (Figure 6). 
Near-isogenic lines carrying the ovate mutation show that shape 
is already determined at anthesis (Van Der Knaap and Tanksley, 
2001) (Figure 6A) and obovoid shape gradually decreases during 
the development of the fruit (Figures 6B,C). 
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The molecular function of OVATE and its family mem- 
bers are not well understood. Yeast two Hybrid (Y2H) screens 
using Arabidopsis KNOX and BELL transcription factors as bait 
led to the identification of OFP members, lending support for 
the notion that OVATE interacts with patterning genes that 
impact fruit shape at the early stages of gynoecium develop- 
ment (Hackbusch et al, 2005; Wang et al., 2010). OFP members 
have also been shown to repress transcription (Wang et al, 2007, 
2011) and overexpression of AtOFPl leads to dwarf phenotypes 
in Arabidopsis and tobacco, in part by negatively regulating the 
transcription of GA20oxl, a key gene in the gibberellin biosynthe- 
sis pathway (Hackbusch et al., 2005; Wang et al., 2007). Contrary 
to findings in Arabidopsis, Y2H of the tomato OVATE protein 
as bait did not lead to the identification of transcription factors 
including KNOX or BELL. Instead, 1 1 out of 26 members of the 
TONNEAU1 Recruiting Motif (TRM) superfamily were identi- 
fied including the putative ortholog of AtTRM 17/20 (Figure 7 
and Table 3). Of all interacting clones obtained, 63.8% belonged 
to the TRM family. The TRM clones identified from the screen 
were partial clones and the overlap between interacting clones of 
the same gene is highlighted in orange (Figure 7). TRMs interact 
with TONNEAUla (TONla), TON lb and TON2/FASS proteins, 
which play critical roles in preprophase band formation and 
microtubule array organization (Camilleri et al, 2002; Azimzadeh 
et al, 2008; Spinner et al, 2010, 2013; Drevensek et al, 2012). 
This finding suggests that OFPs interact with TRMs and micro- 
tubules in addition to acting as transcriptional repressors, and 
thus could provide a mechanistic link between organ pattern- 
ing and growth. TONla, TON lb and TON2 interact with the 
TRM via the M2 and M3 motifs, respectively whereas the TRM 
motif that recognizes OVATE has not yet been identified. Most 
single knockouts of Arabidopsis OFPs exhibit no or mild pheno- 
types (Pagnussat et al, 2007; Li et al., 2011; Wang et al, 2011). 
On the contrary, the premature stop codon mutation found in 
tomato OVATE causes a dramatic morphological change in ovary 
shape, suggesting it may be a unique member of the family. These 
findings together suggest that OVATE acts early in carpel develop- 
ment, possibly during phase 3 corresponding to the link between 



primordia initiation and positioning to growth of the developing 
carpels. 

SUN 

SUN controls fruit elongation, including those found in com- 
mercially grown plum tomatoes, the very long and tapered 
shaped heirloom and oxheart tomatoes (Rodriguez et al., 2011) 
(Figure IB). SUN's effect on fruit elongation is much more pro- 
nounced than the effect of OVATE (Figures 1A,B, 6). The locus 
was fine mapped to the short arm of chromosome 7 and found 
to encode a member of the IQD family of calmodulin-binding 
proteins (Van Der Knaap et al., 2004; Xiao et al, 2008). The muta- 
tion arose from a highly unusual 24.7 kb duplication event from 
chromosome 10 to chromosome 7 (Jiang et al, 2009) (Figure 3, 
Table 2). This transposition was mediated by the retrotranspo- 
son Rider, which has also been found to underlie mutations at a 
few other loci in cultivated tomato unrelated to fruit shape (Jiang 
et al., 2012). Expression of wild type SUN is found in 10 days post 
anthesis fruit but in general is extremely low in all tissues exam- 
ined (Figure 4). The duplication placed SUN in a new genome 
environment leading to much higher expression throughout flo- 
ral and fruit development and an extremely elongated fruit (Xiao 
etal, 2008, 2009). 

The effect of SUN on fruit shape is noticeable at anthesis 
albeit that the effect of the gene is more pronounced immedi- 
ately following fertilization (Van Der Knaap and Tanksley, 2001; 
Xiao et al, 2009; Wu et al., 2011) (Figures 6A,B). The results 
suggest that SUN sets up the patterning before anthesis during 
gynoecium development whereas the execution of the patterning 
plan occurs in part after fertilization. Interestingly, SUN also con- 
trols sepal and terminal leaflet shape and high expression leads 
to twisted stems and leaf rachises (Wu et al, 2011) implying a 
role for this gene in lateral (leaf and sepal) as well as termi- 
nal (fruit) organ development. Epistatic interaction of SUN and 



Table 3 | Tomato TONNEAU1 Recruiting Motif proteins (TRM) that 
interact with OVATE in the Y2H screen. 



Tomato Gene ID 


Arabidopsis 


PBS 2 


Number of 


Percentage 




Ortholog 1 




clones 


of total 


Solyc07g008670.2.1 


TRM5 


A 


31 


16.8 


Solyc09g005750.2.1 


TRM19 


A 


27 


14.6 


Solyc06g083660.2.1 


TRM 17/20 


A 


16 


8.7 


Solyc03g115000.2.1 


TRM3/4 


A 


8 


4.3 


Solyc02g082680.2.1 


TRM26 


A 


8 


4.3 


Solyc09g063080.1.1 


TRM 17/20 


B 


7 


3.8 


SolycOl g094640.2.1 


TRM13/14/15/33 


B 


7 


3.8 


Solyc07g032710.2.1 


TRM30/34 


B 


5 


2.7 


Solyc03g006840.2.1 


TRM25 


C 


6 


3.2 


Solyc08g081160.2.1 


TRM13/14/15/33 


C 


2 


1.1 


Solyc12g007140.1.1 


TRM30/34 


D 


1 


0.5 



' The most likely ortholog(s) in Arabidopsis were determined based on BLAST 
search against TAIR10 Arabidopsis proteins. 

2 PBS, Predicted Biological Score, which is computed to assess the reliability of 
the interaction. A denotes strong and reliable interaction and D denotes weak 
and/or questionable interaction. 
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DNA BD 
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FIGURE 7 | Protein structure of two TONNEAU1 Recruiting Motif (TRM) 
proteins that interact with tomato OVATE in a Y2H study. Full length 
OVATE was used as bait and a truncated form of tomato TRM proteins 
were used as prey. The clones were grown on medium selecting for the 
bait and prey plasmids whereas the X-a-Gal staining highlights the strength 
of the interaction. The most likely Arabidopsis ortholog is listed below the 
gene annotation number. Solyc06g083360 interacts strongly and 
SolycOl g094640 interacts weakly with OVATE. The colored boxes 
designate the M1-M6 motifs defining the TRM family (Drevensek et al., 
2012). The orange domain in the protein is the overlapping region of 
different prey fragments identified in the Y2H screen and the numbers 
indicate amino acid positions. 
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OVATE is likely with respect to growth of the proximal part of the 
fruit (Figure 6C). The degree of obovoid (pear) shape is much 
pronounced in the double NIL than in the sum of the single NILs. 

SUN changes fruit shape by redistributing fruit mass; an 
increase in cells in the proximal-distal direction is accompanied 
by a decrease in cell number in the columella and septum in the 
medio-lateral direction throughout the entire fruit (Wu et al., 
2011) (Figures 6A-C). This suggests that alterations in cell divi- 
sion patterns are critical for fruit shape changes mediated by 
SUN. Yet, how SUN accomplishes changes in cell division pat- 
terns is poorly understood. The IQD members share a common 
central motif of 67 conserved residues named the IQ67 domain 
that binds calmodulin (CaM) (Abel et al, 2005; Levy et al, 2005; 
Huang et al., 2013). High expression of the first identified mem- 
ber of the family, AtlQDl, leads to increases in glucosinolates 
(Levy et al., 2005), a class of secondary metabolites involved 
in plant defense that is absent from Solanaceous plants. How 
increases in glucosinolate levels in Arabidopsis relate to fruit 
shape changes in tomato is therefore, not clear. High expression of 
SUN leads to phenotypes associated with auxin homeostasis, yet 
direct links with auxin through signaling and hormone levels have 
not been established (Wu et al., 2011; Clevenger, 2012). A recent 
Y2H study demonstrated that Arabidopsis IQD1 interacts with 
CaM/CMLs and kinesin light chain-related protein- 1 (KLCR1), 
the latter acts as a motor for transport of vesicles, organelles, 
mRNA-protein complexes within the cytoplasm along micro- 
tubules (Burstenbinder et al., 2013). The directional transport 
of cargo by kinesins could involve the regulation of cell divi- 
sion patterns (Hirokawa et al., 2009; Akhmanova and Hammer, 
2010; Verhey et al, 201 1). The association of AtlQDl with micro- 
tubules suggests that it acts as a scaffold protein to recruit cargo 
to kinesin motors for directional transport along microtubules 
(Burstenbinder et al., 2013). Whether SUN plays a similar role 
as AtlQDl by interaction with KLCR1 proteins is unknown. 
However, the possible involvement in transport of cargo and the 
regulation of cell division patterns would suggest that the mutant 
version of SUN that is highly expressed in developing flowers may 
act as early as stage 3 in organ development, similarly to OVATE 
(Figures 2B,C, Table 1). 

CNR/FW2.2 

The first fruit weight QTL that was cloned from vegetables and 
fruit crops was FW2.2 (Frary et al., 2000). The locus was fine 
mapped to the bottom of chromosome 2 and found to encode 
a member of a novel family of cysteine-rich proteins that share 
the PLAC8 motif (Guo et al., 2010). The family is known to reg- 
ulate cell number, hence the new name for FW2.2-like genes: Cell 
Number Regulator (CNR) (Guo et al., 2010; Guo and Simmons, 
2011) (Figure 3). The underlying mutation to cause changes in 
fruit weight was predicted to be in the promoter as there were 
no polymorphisms in the coding region of the gene (Frary et al., 
2000). Association mapping led to the identification of a puta- 
tive promoter mutation that underlies the fruit weight changes 
(Figure 3). Expression of CNR/FW2.2 is in general low, except 
in the root, young flower buds and developing fruit (Figure 4A). 
Its expression is also the highest in the IM/FM reducing to 
nearly undetectable levels in the floral buds 6 days after initiation 



(Figure 4B). The allele increasing fruit weight causes the enlarge- 
ment of the placenta and columella regions of the fruit (Cong 
et al., 2002; Gonzalo et al., 2009). Previous studies suggested that 
the members of the CNR family are localized to the membrane 
facilitating the transport of ions such as cadmium (Song et al., 
2004) and calcium (Nakagawa et al., 2007) across membranes 
(Guo et al., 2010; Libault et al, 2010). Very little additional infor- 
mation is known about the function of CNR/FW2.2 and how 
regulation of ion transport would lead to changes in cell division. 
Ovary size is different at anthesis, implying that CNR/FW2.2 acts 
early during development of the gynoecium. Based on expression 
profile, the promoter mutation may result in fruit weight changes 
as early as phase 1 or 2 (Table 1). 

TOMATO FRUIT WEIGHT AND SHAPE GENES ACTING 

P0ST-ANTHESIS 

SUN 

SUN clearly impacts the patterning of the fruit prior to anthesis 
(see above). However, the most dramatic effect of SUN on shape 
is manifested after anthesis, during phase 5, which is the cell divi- 
sion stage of fruit development (Van Der Knaap and Tanksley, 
2001; Xiao et al, 2009) (Figure 6B). As a result of SUN expres- 
sion, cell number was much higher along the proximal-distal axis 
and lower along the medio-lateral axis at 7 days post anthesis 
compared to anthesis (Wu et al., 2011) which are likely due to 
the changes in cell division rates in one direction over another 
and not the duration of cell division since fruit ripening time 
is not altered (data not shown). The proposed changes in cell 
division rates in different tissues of the developing fruit is likely 
because fruit weight is not altered and thus SUN appears to result 
in a redistribution of mass. This change in shape is accompa- 
nied by changes in gene expression profiles that are specific to the 
developing pericarp and columella, especially for genes related to 
cell division (Clevenger, 2012). These findings suggest that the 
differences in growth along the various axes after anthesis are 
accompanied by differential gene expression to achieve the final 
fruit shape. These differences in gene expression in the different 
tissue types cede at the time when fruit shape mediated by SUN is 
final which is around 10 days post-anthesis (Clevenger, 2012). 

SLKLUH/FW3.2 

The second fruit weight QTL identified from vegetable and fruit 
crops is FW3.2 (Chakrabarti et al., 2013). The gene was fine 
mapped to the bottom of chromosome 3 encoding a cytochrome 
P450 of the CYP78A class and the likely ortholog of Arabidopsis 
KLUH (Zhang et al., 2012; Chakrabarti et al., 2013) (Figure 3). 
Based on association mapping and additional segregation exper- 
iments, a mutation in the promoter of SIKLUH is proposed 
to underlie the change in tomato fruit weight. This mutation 
is located 512 bp upstream of the predicted start of SIKLUH 
transcription in a putative cz's-element that is known as an organ- 
specific element found in nodulin and leghemoglobin genes 
(Stougaard et al., 1990; Chakrabarti et al, 2013) (Figure 3). 
Expression of tomato KLUH is high in young growing tissues 
containing meristems or developing seeds (Figure 4A). Also, its 
expression is particularly high in the IM/FM and decreases in the 
developing flower buds (Figure 4B). Moreover, within the fruit, 
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KLUH is very highly expressed in the developing seeds and much 
lower in the developing pericarp (Chakrabarti et al., 2013). 

The mutant allele of SIKLUH, found in many cultivated tomato 
accessions, does not impact ovary size at anthesis; rather its effect 
on fruit weight becomes apparent 3 weeks post-anthesis (Zhang, 
2012). Yet, transgenic down regulation of SIKLUH led to shorter 
plants and leaves, smaller flowers in addition to reduced fruit 
weight (Chakrabarti et al., 2013). This result implies that the role 
of KLUH in plant development is broader than the differences 
in the function of the natural KLUH alleles demonstrate. The 
increase in fruit weight arises primarily from increased pericarp 
and septum areas, resulting from additional number of cells. The 
increases in cell number is likely the result of a change in duration 
of cell division and not the rate since fruit ripening is delayed as 
well. In addition to fruit weight, SIKLUH has a pleiotropic effect 
on branching behavior. The large fruit allele of SIKLUH causes 
reduced branch number and length as well as fewer fruits. This 
leads to comparable yields from NIL plants carrying the wild type 
or the mutant SIKLUH allele (Chakrabarti et al, 2013). 

It has been hypothesized that KLUH generates a mobile growth 
promoting signal different from the known phytohormones. 
However, the exact molecular and biochemical nature of the 
"mobile" signal remains elusive and the substrate for this subfam- 
ily of P450 enzymes is also yet to be deciphered (Anastasiou et al., 
2007; Adamski et al, 2009). 

DO 0RTH0L0GS OF TOMATO FRUIT WEIGHT AND SHAPE 
GENES IMPACT FRUIT MORPHOLOGY IN OTHER 
DOMESTICATED PLANTS? 

The domestication of fruit and vegetable crops was likely driven 
by selections for increases in fruit weight and shape in many 
incipient crop species. Thus, the question arises whether any of 
the tomato genes or members of their families are associated with 
fruit weight and shape in other species. Of the fruit weight genes, 
other members of the CYP78A class to which SIKLUH/FW3.2 
belongs are known to regulate floral organ and fruit size, leaf 
and seed size, embryo and endosperm size, apical dominance 
and plastochron length in Arabidopsis, moss and rice (Ito and 
Meyerowitz, 2000; Miyoshi et al, 2004; Anastasiou et al, 2007; 
Adamski et al, 2009; Katsumata et al, 2011; Fang et al, 2012; 
Nagasawa et al., 2013). More intriguingly, in Capsicum spp (chile 
pepper), Cucumis melo (melon) and Vitis vinifera (grape), the 
putative ortholog of KLUH and members of the same CYP78A 
class were associated with larger fruit, suggesting a possible role 
of this small and largely unknown cytochrome P450 family in 
parallel domestication processes in fruit and vegetable crops 
(Chakrabarti et al, 2013; Doligez et al., 2013; Monforte et al, 
2014). Collectively, these findings point toward an evolutionarily 
highly conserved function for this subfamily of P450s in regu- 
lating plant organ size. For CNR/FW2.2, members of the family 
regulate plant growth and biomass as well as ear length and ker- 
nel number per row in maize (Guo et al., 2010) and the number 
of nitrogen-fixing nodules in soybean (Libault et al., 2010). QTL 
studies into the regulation of fruit weight in chile pepper, melon 
and cherry have also implied a possible role for FW2.2/CNR-like 
genes to control weight in a range of crop species (Paran and 
Van Der Knaap, 2007; De Franceschi et al., 2013; Monforte et al., 
2014). 



Of the fruit elongation genes, down regulation of a mem- 
ber of the OFP family in pepper led to a longer shaped fruit 
(Tsaballa et al., 2011), whereas in melon several OFP members 
mapped to fruit shape QTLs (Monforte et al., 2014). This suggests 
that the OFP family is likely to control shape of other fruit and 
vegetables. Of the locule number genes, a weakly overexpressed 
WUSCHEL-like gene in soybean showed an enlarged gynoecium 
(Wong et al, 2011) which also implies that natural alleles of 
WUS could impact the size of fruits and vegetables in other 
crops. 

CONCLUSIONS 

Recent discoveries have started to shed light on the regulation of 
fruit shape and weight, and the molecular mechanisms underly- 
ing this diversity found in cultivated germplasms. However, these 
six genes are unlikely to represent the entire repertoire of genes 
acted on by domestication and diversification. The identification 
of suppressors of ovate (Rodriguez et al, 2013) and the effects 
of genetic background on the severity of the Ic and fas mutants 
both provide evidence for the existence of other genes that inter- 
act with these major regulators of fruit shape and size. In addition, 
the identification of additional fruit weight QTLs (Huang and Van 
Der Knaap, 2011) will result in the identification of new regula- 
tors in fruit weight. Further, the exploitation of TILLING mutants 
that impact shape and weight may also significantly augment the 
resources available in the fruit morphology tool kit (Okabe et al, 
2011). The molecular and biochemical characterization of the 
genes and encoded proteins in the future will greatly add to our 
understanding into the pathways regulating the final dimensions 
of the fruit. 

Advancing the research into the function of fruit morphology 
proteins is going to lead to fundamental insights into plant devel- 
opmental processes. Especially processes that regulate cell prolif- 
eration and enlargement patterns, as well as its rate and duration 
are of particular importance since they pertain to growth of all 
plant organs and eventually yield. In all, the discoveries made 
using tomato fruit morphology as a model will undoubtedly sup- 
port fundamental and applied research that is applicable to many 
other plant systems. 
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